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EXECUTIVE  SUMMARY 


This  report  summarizes  characteristics  of,  and  some  preliminary  results 
from,  a model  developed  to  simulate  fluctuations  In  acoustic  transmission  Induced 
by  Internal  wave  activity  In  the  ocean.  An  explicit  objective  of  this  work  was 
to  provide  a mechanism  for  exploring  transmission  fluctuations  which  might  affect 
the  performance  of  low  frequency  sonar  systems.  The  principal  output  of  the 
model  Is  simulated  sample  functions  of  transmission  parameters  (amplitude, 
phase)  versus  time,  for  specified  frequencies  and  source  receiver  geometries. 

These  time  series  can  either  be  directly  played  through  a particular  sonar 
processing  algorithm  to  provide  simulated  sonar  behavior,  or  appropriate  statistical 
descriptors  can  be  extracted  from  the  series,  and  used  to  analytically  estimate 
sonar  performance. 

The  Internal  wave  field  description  In  the  model  Is  based  on  a power 
spectrum  proposed  by  Garrett  and  Munk.  This  spectrum  Is  used  to  generate  explicit 
sample  functions  of  Internal  wave  activity,  evaluated  over  a selected  mesh  of 
space  time  points.  From  this  set  of  values,  Internal  wave  Induced  perturbations 
In  the  speed  of  sound  at  these  points  may  be  computed.  The  effect  of  these  per- 
turbations on  acoustic  propagation  Is  then  estimated  through  a ray  tracing 
approximation.  The  basic  assumption  Is  that  the  Internal  wave  effects  cause 
significant  variations  only  In  the  phase  of  ray  paths  connecting  source  and 
receiver.  As  the  Internal  wave  field  evolves  In  time,  these  phases  undergo 
complex  fluctuations  whose  nature  depends  on  the  particular  geometry  of  each 
ray  path.  In  turn,  the  coherent  sum  of  these  ray  paths  (e.g.  the  sound  energy 
observed  at  a receiver)  fluctuates  with  time. 

While  the  basic  form  of  the  Internal  wave  spectrum  used  Is  general, 
some  of  the  parameters  Involved  depend  on  characteristics  of  the  local  sound 
speed  profile.  In  addition  the  ray  path  parameters  used  In  the  model  depend 
In  detail  on  this  profile  and  other  environmental  parameters.  Hence  the  model 
Is  suitable  for  Investigating  the  sensitivity  of  computed  fluctuations  to 
differing  environmental  conditions.  Also,  the  current  Implementation  differs 
from  some  previous  models  In  that  It  Is  not  restricted  to  cases  where  the 
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direction  of  Internal  wave  propagation  parallels  that  of  acoustic  propagation; 
Instead  the  more  reasonable  assumption  of  Isotropy  In  propagation  direction 
Is  adapted  for  the  Internal  waves.  Finally,  as  already  noted,  since  the  model 
simulates  actual  time  series,  rather  than  calculating  fixed  analytic  parameters 
of  such  series  directly.  It  offers  considerable  flexibility  In  applications  to 
sonar  system  assessment. 

Preliminary  runs  of  the  model  have  been  made  for  four  different  sets 
of  environmental  conditions.  Of  these,  the  first  two  sets  were  designed  to 
explore  the  sensitivity  of  model  results  to  changes  In  Input  parameters.  The 
first  of  these  environments  Is  nominally  typical  of  the  Atlantic  Ocean  near 
Bermuda  In  autumn,  while  the  other  represents  the  North  Atlantic  In  summer. 

The  other  two  sets  of  runs  were  explicitly  selected  to  match  conditions  observed 
dur»ng  two  measurements  of  CW  transmission  fluctuations:  one  an  experiment  done 
near  Eleuthera,  and  the  other  an  experiment  conducted  over  an  Eleuthera  Bermuda 
Transmission  Path. 

In  the  first  two  cases  runs  were  made  at  frequencies  of  50  and  220  Hz 
for  source-receiver  separations  ranging  from  5 to  200  nml , and  for  several 
source-receiver  depth  combinations.  Figure  1 shows  a typical  time  series 
resulting  from  these  runs.  Figure  2 shows  plots  of  statistical  parameters 
(standard  deviation  and  relaxation  time)  versus  source-receiver  separation 
for  transmission  loss  time  series  In  the  Bermuda  (autumn)  environment,  at 
220  Hz,  and  at  50  Hz.  The  large  scatter  In  the  50  Hz  results  appears  to  be  a 
consequence  of  statistical  uncertainties,  but  the  220  Hz  results  are  consistent 
with  a hypothesis  that  the  standard  deviation  should  Increase  with  range  to  a 
"saturation"  value  of  5 to  6 dB,  and  then  remain  constant,  while  the  relaxation 
time  should  decrease  with  range. 

Figure  3 shows  a comparison  of  406  Hz  propagation  fluctuation  measure- 
ments made  near  Eleuthera,  with  statistics  derived  from  a model  run  made  using 
environmental  parameters  designed  to  match  the  experimental  conditions.  The 
figure  plots  the  rate  (crossings  per  hour)  at  which  the  time  series  crosses  a 
specified  level  as  a function  of  that  level.  The  horizontal  extent  of  this 
plot  Is  a measure  of  the  magnitude  of  fluctuations,  while  the  height  reflects 
the  fluctuation  rates  Involved.  The  agreement  between  data  and  model  Is  considered 
very  good. 

More  generally  the  model  results  match  typical  experimental  observations 
of  CW  transmission.  Of  these,  the  most  striking  are  frequent  deep  fades  In 
propagation  amplitude,  and  the  character  of  related  phase  fluctuations  which 
take  place  relatively  slowly,  but  with  large  (many  cycles)  amplitude. 

The  ensemble  of  model  runs  made  to  date  Indicate  a number  of 
trends  which  seem  to  be  significant.  These  may  be  summarized  as  follows: 

e The  rate  at  which  propagation  fluctuations  occur,  for 
a fixed  source  and  receiver,  Increases  (the  relaxation 
time  decreases)  as  the  acoustic  frequency  Increases 

e There  Is  some  evidence  that  fluctuation  rates  Increase 
with  Increasing  source-receiver  separation 
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FIGURE  3 AXIS  CROSSING  RATE  FOR  ElEUTHERA  DATA  AND  MODEL  PREDICTIONS 
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• Both  predicted  and  observed  fluctuations  tend  to  a 
maximum  value  of  about  5 to  6 dB,  a value  which  Is 
characteristic  of  phase  random  multi-path  propagation. 

One  set  of  runs  shows  evidence  that  this  value  decreases 
with  range  at  short  (less  than  40  nml)  separations  but 
the  current  runs  are  Insufficient  to  determine  whether 
this  feature  changes  with  source-receiver  depth,  frequency, 
or  environment 

• Computed  relaxation  times  range  from  one  to  two  minutes 
up  to  about  30  minutes.  Values  at  220  Hz,  at  ranges 

of  100  to  200  nml,  are  about  five  minutes  In  both  the 
Bermuda  and  Atlantic  envl ronments , and  on  the  order  of 
one  to  two  minutes  In  the  two  400  Hz  cases  examined. 

Conventionally,  sonar  fluctuation  parameters  are  taken  to  be 
Independent  of  range  and  essentially  Independent  of  frequency  over  the  operating 
range  of  any  particular  sonar.  Empirical  fits  to  low  frequency,  long  range, 
sonar  performance  suggest  relaxation  times,  for  the  detecting  process,  of  hours 
rather  than  minutes.  The  possibility,  suggested  by  the  results  cited  In  this 
study*  that  fluctuation  parameters  may  vary  with  both  frequency  and  range  could 
be  of  considerable  Importance  In  optimizing  sonar  design  and  use.  Crudely 
speaking,  the  cited  propagation  fluctuation  characteristics  Imply  relatively 
better  performance  with  Increasing  range  and  frequency,  than  Is  predicted  by 
conventional  calculations.  For  example,  to  the  extent  that  such  effects 
dominate,  sonars  designed  to  reach  a given  range  could  have  lower  gains,  and 
operate  at  higher  frequencies,  then  would  otherwise  be  necessary. 

Of  course,  many  phenomena  besides  those  treated  can  produce  or 
modify  observed  fluctuations  on  any  particular  sonar.  In  a complete  description 
these  must  be  Included  also.  In  addition  the  sonar  system  Itself  modifies  the 
properties  of  an  Incoming  signal  so  that,  ultimately,  the  observed  fluctuations 
depend  on  the  properties  of  the  system.  Long  Integration  times  will  remove,  or 
smooth,  relatively  rapid  fluctuations.  High  resolution  systems  can  Introduce 
fluctuations  by  responding  to  changes  In  parameters  which  would  not  be  resolvable 
on  less  powerful  systems.  This  Is  not  to  suggest  that  the  trends  observed  In 
this  study  are  not  of  Interest,  but  only  that  realistic  assessment  of  sonar 
performance  must  Integrate  such  results  Into  a more  comprehensive  sonar  model. 

As  such  models  are  developed,  estimates  of  transmission  statistics  will  play 
a central  role. 
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I.  INTRODUCTION 


Modeling  and  assessment  of  the  performance  of  narrowband  sonar 
systems  requires  information  on  both  the  mean  value,  and  the  variability, 
of  acoustic  Irputs  (signal  and  noise)  to  the  system.  While  substantial 
attention  has  been  given  to  the  problem  of  estimating  the  mean  values  from 
environmental  and  system  parameters,  the  problem  of  predicting  corresponding 
statistical  properties  has  been  given  far  less  emphasis.  This  report  con- 
siders, in  particular,  the  question  of  time  variations  in  transmission  loss 
between  a fixed  source,  and  a fixed  receiver  in  the  ocean.  There  are  two 
bases  for  this  choice  of  topics.  First,  transmission  loss  is  a critical 
element  in  determing  both  received  signal  levels  and  received  noise  levels. 

In  addition,  there  now  exists  a considerable  body  of  information  on  details 
of  acoustic  transmission  phenomena,  some  of  which  provide  a basis  for  pre- 
dicting stochastic  elements  of  propagation  loss. 

Of  these  phenomena,  this  report  concentrates  especially  on  trans- 
mission loss  perturbations  induced  by  internal  wave  activity.  It  is  generally 
conceded  that  this  is  one  of  the  most  significant  mechanisms  affecting  propa- 
gation at  acoustic  frequencies  between  a few  Hz  and  a few  hundred  Hz.  The 
bulk  of  the  material  which  follows  concerns  development,  and  ini ttal  appl ication 
of  a model  for  simulating  acoustic  propagation  through  an  internal  wave  field. 

The  Internal  wave  description  used  rests  on  a power  spectrum  proposed  by 
Garrett  and  Munk,  and  assumes  a horizontally  isotropic  distribution  of  internal 
waves.  The  acoustic  propagation  calculations  use  ray  tracing  techniques,  and 
require  propagation  parameters  (velocity  profile,  bottom  class,  etc.)  appropriate 
to  the  particular  problem  of  interest.  The  principal  output  of  the  model  is 
estimated  transmission  loss  versus  time.  Such  simulated  time  series  may  be 
used  directly  (through  additional  simulations),  or  indirectly  (through  cal- 
culated values  of  statistical  parameters),  to  evaluate  the  expected  impact 
of  predicted  fluctuations  on  sonar  performance. 
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In  the  next  section,  we  review  briefly  the  elements  of  the  transmission 
i fluctuation  program.  Including  discussion  of  some  of  the  Important  perturbing 

S mechanisms.  In  Section  III,  we  discuss.  In  qualitative  terms,  the  nature  of 

Internal  waves,  and  how  we  model  their  effect  on  acoustic  transmission.  (A 
detailed  quantitative  discussion  of  these  Issues  Is  given  In  the  Appendices.) 

? Section  IV  displays  and  discusses  results  from  the  model  developed  In  this 

* study,  and  Section  V summarizes  study  conclusions  and  recommendations. 
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BACKGROUND 


2.1  ACOUSTIC  FLUCTUATIONS  AND  SONAR  PERFORMANCE 

Observation  of  the  outputs  of  low  frequency  passive  sonars  readily 
establishes  that  the  signals  In  such  systems  normally  exhibit  very  large 
random  fluctuations.  Such  fluctuations  arise  partly  because  the  output 
levels  of  sound  sources  In  the  ocean  vary  with  time,  partly  because  of  the 
way  sound  propagates  from  the  sources  to  the  sonar  receiver,  and  partly 
because  of  attributes  of  the  sonar  system  Itself.  In  particular,  the 
propagation  characteristics  of  the  medium  contribute  to  fluctuations  In 
two  ways:  first,  when  either  source  or  receiver  are  In  motion,  spatial 
variations  In  transmission  properties  will  cause  temporal  variations  In 
received  acoustic  energy;  second,  temporal  variations  In  transmission 
properties  will  cause  fluctuations  In  the  observed  acoustic  fields,  even 
when  source  and  receiver  are  stationary.  In  this  report,  we  focus  par- 
ticularly on  this  last  effect. 

While  the  problem  of  estimating  the  mean  values  of  acoustic  signals 
received  by  sonar  systems  has  received  continuing  and  aggressive  attention, 
the  problem  of  predicting  the  statistical  properties  of  these  signals  has 
been  given  far  less  emphasis.  In  1975  Code  431  of  the  Office  of  Naval 
Research  and  the  Naval  Ship  Research  and  Development  Center  jointly  sponsored 
a symposium  on  fluctuations  and  sonar  performance  which  explicitly  addressed 
how  temporal  fluctuations  affect  passive  sonar  system  performance.1  This 
symposium  reviewed  a considerable  range  of  topics  dealing  with  signal  and 
noise  temporal  fluctuations,  and  the  formalisms  available  for  Incorporating 
such  phenomena  Into  passive  sonar  system  performance  predictions.  The  prin- 
cipal conclusion  of  the  symposium  was  that,  although  fluctuation  phenomena 
are  frequently  a critical  element  in  determining  sonar  performance,  des- 
criptions of  the  significant  phenomena,  and  related  sonar  performance  pre- 
diction methodologies, are  frequently  Inadequate  and  occasionally  nonexistent. 


The  potential  Importance  of  properly  treating  fluctuations  may  be 
readily  appreciated  by  considering  a simple  case:  a narrowband  passive  sonar 
presented  with  fixed  mean  values  of  signal  and  noise  at  Its  Inputs.  In 
such  a system  the  signal  detection  criterion  Is  usually  equivalent  to  asking 
whether  the  signal  to  noise  ratio  exceeds  a specified  threshold.  For  any 
given  mean  value  the  probability  of  this  occurring, the  time  over  which  the 
condition  will  last  (holding  time),  and  the  time  between  occurrences  (recontact 
time),  are  all  determined  solely  by  fluctuations  In  the  Input  quantities. 

When  deterministic  variations  (e.g.,  due  to  target  motion  or  diurnal  changes) 
are  added  the  effect  of  random  fluctuations  cannot  be  separated  out  so 
readily,  but  the  contribution  to  system  performance  Is  no  less  Important. 

The  complexity  of  the  fluctuation  problem  may  also  be  readily 
appreciated.  Sonar  Inputs  are  affected  by  a number  of  processes,  each  with 
Its  own  set  of  stochastic  parameters.  The  resulting  composite  fluctuations 
are  superpositions  of  many  distributions  and  time  scales.  Furthermore,  the 
character  of  the  observed  fluctuation  is  affected  not  least  by  the  properties 
of  the  sonar  system  itself:  bandwidth,  integration  time,  processing  metho- 
dology and  sensor  directivity  all  help  shape  how  the  observed  signals  vary. 
Ultimately,  therefore,  the  character  of  sonar  fluctuations  is  a function  as 
much  of  the  sonar  system  of  interest  as  it  is  of  the  external  acoustic  environ- 
ment. 

2.2  SOUND  SPEED  MICROSTRUCTURE 

In  this  report  we  deal  particularly  with  the  assessment  of  time 
dependence  of  low  frequency  propagation  loss  In  the  ocean.  The  Interest 
in  this  topic  reflects,  of  course,  the  fundamental  Importance  which  proper 
treatment  of  transmission  loss  has  in  determining  sonar  detection  performance. 
In  the  familiar  context  of  the  sonar  equation,  transmission  loss  enters  In 
a number  of  ways:  explicitly  in  assessing  the  loss  from  target  to  receiver, 
and  Implicitly  In  estimating  such  factors  as  ambient  noise,  array  gain, 
reverberation  levels,  etc.  Standard  treatments  of  the  sonar  detection 
problem  treat  all  these  factors  as  Independent,  but  they  are,  in  fact,  linked 
by  common  transmission  loss  phenomena,  and  the  ability  to  model  sonar  detection 
should  be  considerably  Improved  If  the  statistics  of  this  common  factor  can 
be  better  understood. 

In  discussing  acoustic  transmission  phenomena.  It  is  necessary  to 
decide  In  what  terms  the  discussion  will  take  place.  Although. In  principle, 
acoustic  transmission  In  the  ocean  must  be  treated  as  a wave  phenomenon,  it 
Is  well  known  that  open  ocean  conditions  are  very  often  such  as  to  allow 
approximating  transmission  effects  through  various  ray  trace  techniques. 

Such  techniques  permit  estimates  of  the  two  most  Important  characteristics 
of  open  ocean  transmission:  refraction  effects  due  to  a spatially  varying 
Index  of  refraction  (sound  velocity),  and  reflection  and  absorption  at  the 
ocean  surface  and  bottom.  Such  calculations  properly  deal  with  properties 
of  the  medium  whose  spatial  variation  occurs  on  a scale  at  least  several 
times  larger  than  an  acoustic  wave  length  and  lead  to  the  familiar  picture 
that  sound  transmission  between  two  specified  points  takes  place  along  a 
number  of  well  defined,  but  distinct  paths  (rays). 


Of  course,  the  paths  which  result  from  ray  tracing  do  actually  cor- 
respond to  approximate  solutions  of  the  wave  equation.  The  wavefronts 
(regions  of  constant  phase)  for  each  ray  are  the  surfaces  perpendicular 
to  the  ray.  Phase  shifts  along  a given  path,  for  any  frequency,  can  be 
calculated  from  the  sound  speed  at  each  point  along  the  path.  With  this 
Interpretation  many  nominally  wave  effects  can  be  estimated  from  ray  trace 
results.  Interference  effects  are  obtained  by  adding  multipaths  coherently 
using  the  intensities  and  phases  derived  from  ray  tracing.  Responses  of 
extended  sets  of  transducers  (arrays)  may  be  computed  using  phases  similarly 
obtained.  Under  many  circumstances  certain  diffraction  effects  can  be  treated 
as  a straightforward  perturbation  of  the  ray  trace  solutions.  In  a similar 
manner.  It  Is  often  possible  to  Interpret  and  evaluate  those  mechanisms  which 
cause  transmission  loss  fluctuations  in  terms  of  a ray  trace  picture. 

It  is  readily  apparent  that,  from  this  viewpoint,  there  are  two 
classes  of  fluctuation  phenomena:  those  In  which  the  amplitudes  of  the  compo- 
nent multipaths  vary  with  time,  and  those  in  which  their  phases  vary.  While 
the  latter  effect  does  not  change  the  intensities  of  individual  ray  paths, 

It  can  by  altering  their  relative  phases,  cause  variations  In  the  intensity 
of  their  sum,  and  hence  generate  fluctuations  in  overall  transmission  loss. 

The  physical  phenomena  which  cause  these  variations  can  also  be 
divided  into  two  principal  catagorles:  interactions  of  ray  paths  with  the 
time  varying  sea  surface,  and  propagation  of  ray  paths  through  structure 
within  the  ocean  which  changes  with  time.  The  first  of  these  involves 
fluctuations  which  arise  as  ray  paths  reflect  from  the  corrugated  moving 
boundary  presented  by  surface  waves.  The  second  involves  a number  of 
different  phenomena  each  of  which  contribute  to  time  variations  In  the 
ocean's  transmission  properties. 

The  starting  point  for  ray  trace  calculations  is  the  ocean  sound 
speed  profile:  the  functional  dependence  of  sound  speed  in  the  sea  on  depth 
and  position.  Typically  these  profiles  characterize  sound  speed  dependence 
with  spatial  resolutions  of  several  tens  of  meters  in  depth  and  several 
hundreds  to  several  thousands  of  meters  in  range.  Careful  observations 
establish,  however,  that  the  sound  speed  at  any  point  is  not  constant,  and 
fluctuations  occur  on  spatial  scales  which  range  from  much  finer  than,  to 
much  grosser  than,  the  nominal  resolutions  cited  above.  These  fluctuations 
act  to  perturb  acoustic  transmission  paths,  causing  temporal  variations  In 
transmission  loss. 

Much  of  the  observed  fluctuation  in  sound  speed,  in  turn,  arises 
from  mechanisms  which  cause  mixing  of  ocean  waters.  Sound  speed  In  the 
ocean  Is  a function  of  water  temperature.  In  regions  where  there  is  a 
temperature  gradient,  such  mechanisms  mix  waters  of  different  temperatures. 
This  local  intrusion  of  slightly  different  temperatures  manifests  itself 
as  stochastic  variations  in  the  sound  speed  structure,  which  can  affect 
propagation  characteristics. 
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At  the  largest  scale,  ocean  water  flow  Is  In  the  form  of  global 
currents,  large  scale  eddies,  and  tidal  effects.  Such  motion  can  generate 
large  turbulent  vortices  which  dissipate  their  energy  by  generating  somewhat 
smaller  scale  vortices,  starting  a chain  which  ends  with  very  small  scale 
motions.  The  expected  distribution  over  spatial  scales  in  such  situations 
was  derived  by  Kolmogorov2,  who  predicted  that  the  energy  density  In  tur- 
bulent flow  varies  as  wavelength  to  the  five-thirds  power.  This  spectrum 
of  fluid  movement  Induces  a corresponding  spatial  continuum  of  temperature 
fluctuations,  and  hence  of  sound  speed  fluctuations. 

In  a homogeneous  medium  this  turbulent  spectrum  would  dominate. 

The  real  ocean  Is  stratified,  exhibiting  variations  In  density  and  temperature 
with  depth.  This  stratification  allows  other  types  of  fluid  mixing,  which 
modify  the  purely  turbulent  spectrum.  The  most  notable  of  these  are  Internal 
waves,  which  are  the  major  subject  of  this  study.  Internal  wave  effects 
tend  to  dominate  the  low  wave  number  region  of  the  fluid  velocity  spectrum. 
More  generally  Mosely  and  Del  Balzo3  suggest  that  the  spectrum  of  fluid 
velocity  can  be  divided  Into  seven  regions:  beyond  the  very  large  scale 
currents  which  derive  their  energy  from  geophysical  processes,  there  are 
large  stable  eddies,  the  Internal  wave  region,  a buoyancy  region,  a region 
of  inertial  turbulance,  and  finally,  at  very  large  wavenumbers,  a viscous 
dissipation  region,  where  the  fluid  motion  Is  damped,  and  Its  energy  dis- 
sipated. 


For  low  frequency  (less  than  a few  hundred  Hz)  acoustic  propagation 
there  Is  especial  Interest  in  the  Internal  wave  portion  of  this  spectrum. 

First  this  Is  the  most  energetic  region  of  the  spectrum.  In  addition,  the 
acoustic  wave  lengths  at  low  frequencies  are  on  the  order  of  several  meters 
or  more  In  length.  Much  of  the  inhomogeneity  spectrum  (and  especially  the 
upper  region  where  turbulence  dominates)  has  much  shorter  spatial  scales 
than  this.  Under  these  conditions  the  structure  tends  to  act  as  weak 
isotropic  scatterers,  whose  major  effect  is  to  remove  small  amounts  of 
energy  from  a ray  path.  The  internal  wave  region,  on  the  other  hand, 
presents  Inhomogeneities  large  compared  with  the  acoustic  wavelength. 

These  tend  to  act  as  very  forward  scatterers  whose  net  effect  is  to  pre- 
serve ray  paths,  but  alter  their  phase.  This  random  phase  perturbation, 
by  generating  time  varying  interference  effects,  can  have  a strong  effect 
on  total  transmitted  intensity,  when  more  than  one  acoustic  path  contributes. 

There  exists  a considerable  body  of  work  (much  of  it  developed  for 
atmospheric  propagation)  on  propagation  of  energy  through  media  with  random  % 
Inhomogenletles.  Much  of  this  material  (e.g.,  the  classic  work  by  Tatarskii  ) 
assumes  that  the  Inhomogenietles  follow  the  Kolmogrorov  spectrum.  Much 
of  It,  also,  seeks  analytic  solutions  for  particular  statistical  parameters 
under  conditions  (very  large  wave  number,  short  ranges,  homogeneous  ocean) 
which  severely  limit  application  of  results.  In  the  chapter  which  follows 
we  outline  a model  which  assesses  the  effect  of  internal  waves  on  acoustic 
propagation  through  a simulation  approach.  This  has  the  decided  advantage 
of  not  tying  the  model  to  particular  measures  of  acoustic  variability.  In- 
stead, the  model  generates  explicit  sample  time  series,  which  can  be  used 
to  obtain  any  statistical  parameters  of  interest. 
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III.  INTERNAL  WAVE  MODEL-QUALITATIVE  DISCUSSION 


In  tills  chapter  we  give  a qualitative  description  of  a model 
developed  to  allow  simulating  the  effects  of  Internal  wave  fields  on  acoustic 
propagation.  (A  roughly  parallel  treatment  of  the  corresponding  mathematical 
formalism  Is  given  In  Appendix  A.)  Essentially  the  model  uses  estimates  of 
the  space-time  power  spectrum  of  internal  waves  to  build  up  simulated  sample 
functions  of  the  Internal  wave  field  which  can  be  explicitly  evaluated  at 
any  given  space-time  point.  From  this  function,  estimates  of  the  variation 
In  sound  speed,  due  to  Internal  wave  activity,  can  be  formed.  A standard 
ray  trace  transmission  loss  model  Is  then  used  to  establish  expected  trans- 
mission parameters  In  the  absence  of  Internal  waves,  and  also  to  estimate 
how  the  Internal  wave  field  Induced  sound  speed  fluctuations  perturb  these 
parameters.  For  any  one  sample  function  the  perturbations  can  be  evaluated 
as  a function  of  both  time,  and  source-receiver  position.  The  calculation 
can  be  replicated,  to  build  up  statistics,  by  constructing  additional  sample 
functions. 

The  construction  of  Internal  wave  functions  rests  on  a canonical 
description  of  the  ocean's  density-temperature  structure.  There  is  some 
expectation  that  the  form  of  internal  wavefleld  spectrum  is  not  very  sensitive 
to  local  conditions  so  that  use  of  a single  general  description  may  be  satis- 
factory. On  the  other  hand,  the  ray  trace  calculations  and  ev  Nation  of  certain 
Internal  wavefleld  parameters  are  based  on  estimates  of  the  mean  sound  speed  at 
the  site  and  for  the  conditions  of  Interest,  making  the  model  results  particular 
to  those  conditions.  This  makes  it  possible  to  use  the  model  to  explore  whether 
some  sets  of  propagation  conditions  are  more  sensitive  to  internal  wave  effects 
than  others. 

Internal  waves  are  wave  motions  (somewhat  like  surface  waves)  which 
propagate  In  the  body  of  a fluid,  where  there  is  a density  gradient.  Like 
surface  waves,  they  Involve  vertical  motions  of  the  fluid,  which  Induce  local 
mixing  of  the  fluid.  As  noted  In  the  last  chapter,  If  there  Is  a temperature 
gradient  as  well  as  a density  gradient,  the  effect  of  such  mixing  Is  to 
Introduce  local  anomalies  In  the  sound  speed  profile,  which  can  perturb 
acoustic  propagation. 


Recent  work  by  a number  of  workers  has  begun  to  provide  a quantitative 
foundation  for  understanding  and  modeling  these  effects.  In  outlining  this 
progress,  It  Is  convenient  to  divide  the  problem  Into  three  parts:  the  nature 
of  the  Internal  wave  field,  the  effect  of  Internal  waves  on  the  velocity  pro- 
file (Index  of  refraction),  and  the  modeling  of  transmission  fluctuations 
Induced  by  the  velocity  changes. 

Characterizations  of  the  Internal  wave  field  are  based  on  work  by 
Garret  and  Munk5,  which  combines  mathematical  forms  derived  from  solutions 
of  the  wave  equation  with  empirical  fits  to  data.  Briefly  their  description 
has  the  following  properties: 

• Internal  waves  are  a superposition  of  horizontally 
propagating  plane  waves  with  a specified  propagation 
direction,  wave  number,  and  frequency. 

t Internal  wave  propagation  can  be  sustained.  In  two 
dimensions,  for  a considerable  range  of  wave  numbers 
and  frequencies. 

• For  any  particular  density  gradient  there  Is  a set  of 
Internal  wave  solutions  (normal  modes)  at  discrete  (but 
densely  spaced)  combinations  of  frequency  and  wave 
number  which  span  the  range  of  possible  frequency-wave 
number  values. 

• Garrett  and  Munk  derive  equations  for  the  normal  modes 
assuming  a particular  (exponential)  density  function. 

They  note  that  while  the  form  (vertical  dependence)  of 
these  modes  Is  highly  detailed,  an  RMS  average  over  a 
large  number  of  adjacent  modes  (nearly  constant  frequency - 
wave  number)  Is  quite  smooth.  They  postulate  that  the 
averaged  function  Is  also  relatively  Insensitive  to  the 
choice  of  density  gradient. 

• Coupling  among  the  modes  causes  energy  to  be  shared  among 
them.  Garret  and  Munk  derive  a semi -empirical  expression 
for  the  power  spectrum  (e.g.,  the  averaae  energy  density 
per  unit  wave  number-frequency  Interval).  Since  the 
physical  process  Is  stochastic,  the  actual  amplitude  of 
any  mode  at  any  time  Is  a random  variable  whose  mean 
square  value  Is  determined  by  the  power  spectrum.  A 
natural  assumption  Is  that  the  amplitude  Itself  is  Rayleigh 
distributed  with  random  phase  (e.g.,  a Gaussian  Random 
process).  (The  sum  over  closely  spaced  modes  referred 

to  above  has  essentially  the  same  frequency  and  wave 
number  as  the  constituent  modes,  and  Is  again  a Gaussian 
Random  process.) 
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The  vertical  mixing  of  the  water  column  at  any  given  location  Is 
thus  given  by  a superposition  of  a large  number  of  traveling  waves  with  random 
amplitude  and  phase.  The  change  In  temperature  at  that  location  is  assumed 
to  be  proportional  to  the  vertical  displacement  and  the  mean  thermal  gradient, 
and  Is  therefore  Itself  a random  process.  The  change  In  temperature  causes 
a change  In  the  local  speed  of  sound  which  In  turn  Is  equivalent  to  a local 
change  In  the  acoustic  Index  of  refraction. 

In  principle  such  variations  In  the  Index  of  refraction  extend 
vertically  through  the  entire  water  column.  In  fact,  the  Garret  and  Munk 
vertical  distributions  Indicate  that  the  bulk  of  the  effect  takes  place  over 
a relatively  narrow  range  of  depths  at  the  base  of  the  thermocllne.  Porter, 
et  al*,  take  advantage  of  this  by  Introducing  an  approximation  where  the  total 
change  In  Index  of  refraction  (over  all  depths)  Is  imagined  to  be  concentrated 
In  a thin  layer  at  the  base  of  the  thermocllne.  In  this  approximation  the 
variations  In  index  of  refraction  are  still  expressed  as  a superposition  of 
horizontally  propagating  random  waves,  but  the  effect  Is  confined  to  a thin 
slab  In  depth. 

With  this  approximation  It  is  straightforward  to  deal  with  the  effects 
of  Internal  waves  on  acoustic  propagation.  In  the  context  of  ray  trace  treat- 
ments of  propagation  effects.  Any  ray  which  does  not  cross  the  Internal  wave 
layer  Is  unaffected;  a ray  which  does  cross  behaves  like  any  ray  traversing 
a thin  slab  with  Index  of  refraction  different  from  Its  surroundings.  Its 
direction  and  amplitude*  are  essentially  unaffected,  but  its  phase  is  perturbed. 
Thus,  In  this  picture,  the  main  effect  of  internal  waves  Is  to  perturb  the 
phase  of  those  acoustic  paths  which  spend  significant  time  at  the  depths  of 
principal  Internal  wave  activity.  These  perturbations  in  phase  will  cause 
fluctuations  In  acoustic  Intensity  when  ray  paths  are  combined  at  a receiver. 

The  technique  adopted  for  modeling  these  effects  Is  to  use  a Monte 
Carlo  technique  In  which  each  replication  Involves  constructing  a particular 
realization  of  the  internal  wave  field  (e.g.,  picking  amplitudes  and  phase 
for  each  constituent  traveling  wave),  and  then  evaluating  the  resultant 
acoustic  fluctuations  as  a function  of  position  and  time  for  that  particular 
realization.  Data  from  several  such  repl Icatlons  can  be  used  to  provide 
distribution  functions,  and  spatial  and  temporal  correlation  functions .which 
characterize  the  variations  In  acoustic  propagation  which  result  from  the 
Internal  wave  perturbations. 

The  model  developed  here  Is  similar  In  conception  to  models  reported 
on  by  DeFarrarl7  Porter  Splndel  and  Jaffeet  and  Flatte  and  Tappert?  The  work 
by  DeFarrarl  however,  describes  the  Internal  wavefleld  using  only  a single 
mode,  rather  than  the  full  spectrum  Incorporated  In  the  other  models.  The 
work  reported  on  here  closely  parallels  that  of  Porter  et  al.  However,  we 
have  adapted  a considerably  more  sophisticated  acoustic  propagation  module 
than  theirs,  and  we  have  removed  the  simplifying  restriction  that  internal 
waves  propagate  parallel  to  the  acoustic  propagation  direction  only  and 
replaced  It  with  the  more  reasonable  condition  of  isotropy.  Flatte  and 


* There  will  also  be  a change  In  amplitude  due  to  reflections  from  the  layer, 
but  If  the  perturbation  In  Index  of  refraction  is  small,  this  Is  negligible. 
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Tappert  use  parabolic  equation  techniques  to  deal  with  acoustic  propagation 
and  treat  the  Internal  wavefleld  In  three,  rather  than  two,  dimensions.  How- 
ever their  model  deals  with  a canonical  velocity  profile  only,  and  Is  also 
restricted  to  Internal  waves  propagating  parallel  to  the  acoustic  propagation 
direction.  Finally  the  model  developed  here  Is  the  only  one  which  explicitly 
simulates  time  series  of  acoustic  propagation  and  related  parameters.  This 
offers  considerable  power  and  flexibility  In  applications  to  sonar  performance 
problems. 
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IV.  SIMULATION  RESULTS 
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4.1  INTRODUCTION 

In  this  section  we  Investigate  predicted  properties  of  acoustic 
propagation  in  the  presence  of  Internal  waves,  through  use  of  the  simulation 
model  described  In  the  last  section,  and  In  the  appendices  to  this  report. 
Runs  have  been  made  with  this  model  towards  two  distinct  objectives  as 
follows: 

a)  The  acoustic  propagation  elements  of  the  model  deal 
with  specific  sets  of  propagation  parameters  (sound 
speed  profiles,  source/receiver  depths,  water  depth, 
frequency,  etc.).  This  makes  It  possible  to  investi- 
gate how  sensitive  predicted  fluctuations  are  to 
changes  in  these  parameters.  This  has  been  done  for 
two  seperate  sets  of  environmental  parameters  , one 
representative  of  the  Bermuda  area  in  Autumn,  and 
one  representing  the  North  Atlantic  in  Summer.  In 
each  case  runs  were  made  for  several  frequency-source 
receiver  depth-source/receiver  range  combinations. 

b)  Again,  because  of  the  specific  nature  of  the  model,  it 
Is  possible  to  run  direct  comparisons  with  experimental 
data.  Detailed  agreement  in  such  comparisons  Is  hardly 
to  be  expected,  but  gross  features  should  correspond. 

In  any  case,  It  Is  such  comparisons  which  can  most 
effectively  Identify  deficiencies  in  the  model's  predictive 
capabilities.  Below  we  deal  briefly  with  two  such 
comparisons,  one  Involving  data  from  a propagation 
experiment  conducted  off  of  Eleuthera9,  and  another 
from  an  experiment  conducted  between  Eleuthera  and 
Bermuda.10 
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In  each  model  run,  the  basic  result  was  a four  hour  time  series  of 
propagation  loss.  Plots  of  these  allow  qualitative  inspection  of  predicted 
fluctuations.  For  several  runs,  estimates  of  the  time  variation  of  phase 
along  one  or  two  selected  ray  paths  were  also  generated.* 

Finally,  for  each  time  series,  the  mean,  the  variance  and  the 
relaxation  time  of  the  propagation  loss,  in  decibels,  were  estimated.  (Decibel 
values  were  used  since  these  are  normally  the  quantities  of  interest  in  sonar 
performance  calculations.)  The  mean  and  variance  were  computed  using  standard 
definitions.  The  relaxation  time  was  defined  as  the  interval  over  which  the 
autocorrelation  function  for  the  time  series  drops  to  1/e  of  its  initial  value. 
(Note  that  the  mean  of  a set  of  dB  values  is  not  the  same  as  the  mean  of  the 
equivllent  power  values,  even  when  the  latter  is  expressed  in  dB.  Hence  the 
means  quoted  below  will  not  converge  to  conventional  "incoherent"  propagation 
loss  estimates  even  when  predicted  phase  fluctuations  are  large  enough  to 
make  the  constituent  multipaths  essentially  Incoherent.) 

4.2  ENVIRONMENTAL  PARAMETERS 

Figure  4.1  shows  velocity  profiles  for  the  four  sites  of  interest. 

The  first  two  are  from  the  Naval  Oceanographic  Office  "Oceanographic  Atlas,"11 
which  cites  class  averages  for  measured  velocity  profiles  over  selected  regions 
of  the  ocean.  The  last  two  are  profiles  reported  in  the  two  experiments 
cited  above.  Figures  4.2  and  4.3  show  conventional  (Incoherent  sum)  propa- 
gation loss  versus  range  curves  for  the  Bermuda  (autumn)  and  North  Atlantic 
(summer)  environment  as  predicted  by  the  propagation  module  of  the  model, 
in  the  absence  of  internal  wave  activity. 

Table  4.1  shows,  for  each  run  set,  values  used  for  the  parameters 
which  determine  space/time  scales  for  internal  wave  activity.  The  first  two 
columns  are  derived  from  the  sound  speed  profile,  while  the  last  three  are 
obtained  from  the  first  two  (see  Appendix  B).  The  parameter  n0  is  the  peak 
value  of  the  Brunt  Vaisala  frequency  and  defines  the  peak  frequency  for  in- 
ternal wave  activity.  The  parameter  b is  a length  scale  which  defines  how 
rapidly  the  Brunt  Vaisala  frequency  falls  off  with  depth.  The  "In  Layer 
Depth"  In  the  last  column  is  the  depth  at  which  the  model  assumes  internal 
wave  activity  to  be  localized.  Bottom  class  3 was  assumed  for  all  runs  except 
the  North  Atlantic  runs,  which  assumed  bottom  class  4. 

4.3  RESULTS 

4.3.1  Bermuda  (Autumn)  Environment 

Table  4.2  Is  a compilation  of  all  results  for  the  Bermuda  (autumn) 
environment.  The  emphasis  was  on  runs  at  220  Hz,  and  50  Hz,  with  a 100  foot 


* Typically,  the  most  intense  representatives  of  two  distinct  ray  families 
(e.g.,  the  loudest  RSR  and  the  loudest  RR  path)  have  been  selected  for 
presentation.  Those  displayed  are  Illustrative  examples  only;  actual 
model  calculations  allow  up  to  30  paths. 


4-2 


FREQUENCY  - 220  Hz 


FIGURE  4.2  PROPAGATION  LOSS  FOR  BERMUDA  (AUTUMN)  ENVIRONMENT 


TABLE  4.1 

PARAMETERS  FOR  MODEL  RUNS 


CASE 

AXIS  DEPTH  (h) 
(km) 

_1  dc 
c Hi  o 
(Km-1) 

b 

(Km) 

no 

(Rad/Sec) 

IW  LAYER 
DEPTH  (Km) 

Bermuda 

(Autumn) 

1.22 

-.05 

1.45 

.0058 

.72 

N.  Atlantic 
(Summer) 

1.25 

-.043 

2.0 

.0054 

1.0 

"Eleuthera 

Experiment" 

1.30 

-.039 

1.27 

.0052 

.63 

"Eleuthera- 

1.25 

-.029 

1.33 

.0047 

.66 

Bermuda 

Experiment" 
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source,  and  a 2500  foot  receiver.  There  is  a surface  duct,  and  this  choice 
of  depths  forced  propagation  paths  to  go  through  both  the  surface  layer  and 
the  Internal  wave  layer.  The  runs  with  source  and  receiver  at  100  feet  place 
both  in  the  surface  duct,  while  those  runs  with  both  at  2500  feet  define  a 
case  where  both  are  well  below  the  surface  layer  and  the  internal  wave  layer. 
Figure  4.4  plots  the  statistical  parameters  In  Table  4.2  against  range. 

Figures  4.5  through  4.20  show  the  time  series  for  the  220  Hz  runs,  and 
Figures  4.21  through  4.32  show  the  time  series  for  the  50  Hz  runs. 

Qualitatively,  the  time  series  show  features  which  also  characterize 
typical  experimental  results:  frequent  fluctuations  of  several  dB  or  more, 
with  occasional  deep  fades,  where  the  signal  disappears  entirely.  It  Is 
evident  that  the  220  Hz  results  fluctuate  rather  more  rapidly  than  the  50  Hz 
results,  a conclusion  which  might  have  been  anticipated  since  the  size  of 
the  underlying  phase  fluctuations  is  Inversely  proportional  to  acoustic 
wavelength  (Equation  A20). 

Examination  of  Figure  4.4  suggests  that  the  statistical  descriptions 
of  the  220  Hz  results,  at  least,  show  fairly  clear  range  trends.  It  appears 
that  the  standard  deviation  Increases  with  range  until  it  reaches  about  five 
to  six  dB,  after  which  It  holds  constant.  (This  would  accord  well  with 
predictions  that  propagation  fluctuations  should  have  a standard  deviation 
of  5.5  dB  for  long  range  multipath  propagation.12)  The  relaxation  time 
drops  off  smoothly  with  range,  confirming  the  qualitative  impression, 
suggested  by  the  time  series  plots,  that  the  fluctuations  are  more  rapid 
at  longer  ranges.  It  Is  not  clear  whether  the  spread  among  the  three  stan- 
dard deviation  values  (for  different  source-receiver  depth  combinations)  at 
35  miles  Is  statistically  significant;  if  so  it  would  suggest  that  the 
fluctuations  for  a shallow  source  and  receiver  rise  In  amplitude  more  slowly, 
with  range,  than  do  those  In  cases  where  either  source  or  receiver  is  deep. 
Finally,  It  is  worth  noting  that  while  a source  at  35  miles  distance  is  clearly 
within  a convergence  zone  according  to  Figure  4.2,  there  Is  no  significant 
structure  at  this  point  in  the  200  Hz  relaxation  time  vs.  range  curve,  and 
probably  no  significant  departure  from  general  trends  in  the  standard 
deviation  curve. 

In  contrast  to  the  220  Hz  curves,  the  50  Hz  plots  In  Figure  4.4  show 
no  dlscernable  trends  whatsoever.  One  possible  explanation  is  that  the  values 
plotted  are  dominated  by  random  variations— that  Is,  the  time  series  from  which 
these  parameters  were  estimated  were  not  long  enough  to  provide  statistically 
reliable  results.  Figure  4.4  itself  provides  some  internal  evidence  for  this 
possibility;  the  mean  value  of  the  relaxation  times  shown  at  50  Hz  appears  to 
be  three  to  four  times  that  of  the  values  for  220  Hz,  so  that  time  series  of 
a given  length  at  50  Hz  would  have  fewer  degrees  of  freedom,  and  less  statistical 
reliability,  than  similar  series  at  220  Hz.  Direct  examination  of  the  time 
series  themselves  suggests  the  same  conclusion  qualitatively;  major  fluctuations 
rarely  occur  more  than  once  In  a four  hour  period,  and  each  plot  shows  relatively 
few  fluctuation  maxima  of  any  kind,  indicating  that  the  duration  of  the  series 
is  too  short  to  sample  the  full  range  of  fluctuations. 
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FIGURE  4.4  STANDARD  DEVIATION  AND  RELAXATION  TIME  OF  PROPAGATION  LOSS  VERSUS  RANGE 
AT  220  Uz  AND  AT  50  Hz,  FOR  THE  BERMUDA  (AUTUMN)  ENVIRONMENT 


FIGURE  4.5  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  5 NMI 
ANU  A FREQUENCY  OF  220  llz 


FIGURE  4.6  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  lb  i!MI 
AND  A FREQUENCY  OF  220  HZ 


TIME  (HOURS) 

riGURE  4.7  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  2b  NMI 
AIM)  A FREQUENCY  OF  220  IIZ 


BERMUDA  (AUTUMN) 
RANGE  = 35  NMI 
FREQUENCY  = 220  Hz 
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FIGURE  4.9  PROPAGATION  LOSS  TIML  5LRILS  FOR  A RAilGL  OF  3b  rlMI 
AND  A FKLOUEHCY  OF  220  IIZ 


SOURCE  DEPTH  - 100  ft 
RECEIVER  DEPTH  = 100  ft 


FIGURE  4.10  PHASE  FLUCTUAT IOUS  FOR  DOMINANT  EIGLilRAY  AT  A RAilGE  OF 
36  If  MI  AHO  A FREQUENCY  OF  220  IIZ 


MM  I AMD  A FRLQULHCY  OF  UQ  M L 


OF  35  MM  I Afll)  A FliLQULfICY  OF  2Z0  HZ 


BERMUDA  (AUTUMN) 
RANGE  = 100  NMI 
FREQUENCY  = 220  Hz 


Ui 
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FIGURE  4.14  PHASE  FLUCTUATIOilS  DOMIIIA1IT  L IGEIIKAYS  AT  A RAlJGE  UF  1UO  ,IM1 
AM)  A FRElJUEilCY  OF  2ZO  IIZ 


AIU)  A FREQUENCY  OF  220  HZ 


FIGURE  4.  lit  PHASE  FLUCTUATIONS  FOR  DOMINANT  EIGEilRAYS  AT  A RANGE 
OF  100  Hill  AND  A FREQUENCY  OF  220  IIZ 


TIME  (HOURS) 

FIGURE  4.20  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  200  NMI 

AND  A FREQUENCY  OF  200  NZ 


FIGURE  4.21  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  5 IIMI 
AND  A FREQUENCY  OF  SO  HZ 


FIGURE  4.22  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGL  OF  15  NMI 
AilO  A FREQUENCY  OF  50  HZ 


FIGURE  4.23  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  2b  NMI 
AND  A FREQUENCY  OF  50  HZ 


AND  A FKLQUtUCY  OF  50  IIZ 


FIGURE  4.2S  PROPAGATION  LOSS  TIME  SER1LS  FOR  A RANGE  OF  35  NM1 
AIIU  A FREQUENCY  OF  SO  II 1 


AND  A FREQUENCY  OF  DO  HZ 


FIGURE  4.2 0 PHASE  FLUCTUATIONS  FOP  DOMINANT  LIGEilP.AYS  AT  A RA.H.L 
OF  Jb  Hill  AIIO  A FREQUUIC Y OF  bU  HZ 


BERMUDA  (AUTUMN) 
RANGE  = 100  NMI 


FIGURE  4.29  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  100  UNI 

AIIU  A FREQUENCY  OF  SO  HZ 
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FIGURE  4.30  PHASE  FLUCTUATIONS  FOR  DOMINANT  E IGLNRAYS  AT  A RANGL  OF 
100  NMI  AND  A FREQUENCY  OF  SO  HZ 


All)  A FREQUENCY  OF  50  UZ 


FLUCTUATIONS  FOR  DOMINANT  L IGLNRAYS  AT  A RAHGt  OF  100  NM1 
AlIU  A FREQUENCY  OF  SO  NZ 
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However,  direct  examination  of  the  time  series  permits  some  qualitative 
conclusions,  which  are  not  evident  from  Figure  4.4  In  particular,  we  observe 
that 

• As  noted  above,  the  trend  from  50  Hz  to  220  Hz  suggests 
that  the  rate  of  fluctuations  goes  up  (relaxation  time 
goes  down)  as  frequency  Increases 

• At  50  Hz,  as  at  220  Hz,  relaxation  time  decreases  with 
Increasing  range. 


4.3.2  North  Atlantic  (Summer)  Environment 

Table  4.3  summarizes  results  for  the  runs  made  In  the  North  Atlantic 

I (summer)  environment.  Figures  4.33  through  4.38  show  the  time  series  generated 

In  these  runs. 

The  purpose  of  these  runs  was  to  test  whether  predicted  fluctuation 
characteristics  might  vary  significantly  from  one  environment  to  another.  We 
• observe  that  the  220  Hz  parameters  In  Table  4.3  agree  rather  closely  with 

their  equivilents  In  Table  4.2,  with  the  exception  of  the  relaxation  time 

!at  35  nmi , which  is  considerably  shorter  in  the  Atlantic  environments.  Com- 
parison of  the  corresponding  two  time  series  (Figure  4.8  and  Figure  4.33) 
confirms  that  the  Atlantic  case  exhibits  more  rapid  fluctuations.  The  large 
uncertainties  In  these  calculations  make  It  Impossible,  however  to  determine 
whether  this  reflects  a consistent  environmental  difference. 

As  with  the  Bermuda  runs,  there  Is  some  reason  to  question  the  sta- 
tistical reliability  of  the  50  Hz  results,  but  the  tendency  for  lower  frequency 
runs  to  fluctuate  less  rapidlyls  again  dlscernable. 


4.3.3  Experimental  Comparisons 
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Model  runs  have  been  made  for  conditions  which  match  those  of  two 
sets  of  experimental  data.  One  of  these,  reported  by  Porter,  Splndel , and 
Jaffee  we  refer  to  as  "Eleuthera"  data;  the  other  is  deep  ocean  data  from 
Project  Miml  reported  by  Clark  and  Kronengold,  and  referred  to  here  as 
"Eleuthera-Bermuda"  data.  Both  experiments  were  performed  over  deep  ocean 
paths  using  a 406  Hz  tonal  generated  by  a bottom  mounted  projector  at  Eleuthera. 
The  Eleuthera-Bermuda  results  were  recorded  by  fixed  hydrophones  at  Bermuda, 
and  at  an  intermediate  position  between  Eleuthera  and  Bermuda.  The  Eleuthera 
results  were  recorded  by  free  drifting  sonobuoys.  (A  doppler  tracking  system 
was  used  to  follow  the  buoy  motion  so  as  to  allow  corrections  for  position 
changes  to  be  made.) 

The  model  run  for  the  Eleuthera-Bermuda  data  was  designed  for  com- 
parison with  the  experimental  results  from  the  intermediate  hydrophone  at  a 
range  of  120  nmi;  the  time  series  generated  for  this  case  are  shown  in  Figures 
4.39  and  4.40.  The  Eleuthera  model  run  assumed  a range  of  113  miles;  the 
resulting  time  series  for  that  set  of  conditions  are  shown  in  Figures  4.41 
and  4.42. 
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-100 


AilD  A FREQULfO  OF  22 0 ||Z 


nr, URL  4.34  PHASE  FLUCTUATIONS  FOR  UOMIUAUT  LIGLHRAYS  AT  A RAHUL 
OF  35  HMI  A.ID  A FREiJULIICY  OF  220  HZ 


FIGURt  4. 3G  PR0PAGAT10I1  LOSS  TIME  SERIES  FUR  A R.IAGL  OF  IOO  ilMI 

AilO  A FREQUEilCY  OF  220  HZ 


SOURCE  DEPTH  = 100  ft 
RECEIVER  DEPTH  = 4200  ft 


FIGURE  4.37  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  3S  UHI 

AND  A FREQUENCY  OF  SO  117 


FICURC  4.38  PHASt  FLUCTUATIONS  FOR  UOM1NA.IT  LIGLNRAYS  AT  A RANCL 
OF  3b  I1MI  AND  A FRLQUCNCY  OF  bO  IIZ 
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ELEUTHERA 
RANGE  = 113  NMI 
FREQUENCY  = 406  Hz 


FIGURE  4.41  PROPAGATION  LOSS  TIME  SERIES  FOR  A RANGE  OF  113  NMI 
AND  A FREQUENCY  OF  406  HZ 
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Mot  surprisingly,  the  two  time  series  are  very  similar  in  their 
properties,  with  essentially  the  same  mean  loss  (101.7  dB)  nearly  the  same 
a (approx1m»tely  5 dB),  and  nearly  the  same  relaxation  time  (approximately 
une  minute;,  inis  reflects  the  very  similar  conditions  (envi ronment,  range, 
frequency,  source  depth,  receiver  depth)  prevailing  in  the  two  experiments. 

Experimental  results  for  the  Eleuthera-Bermuda  experiment  were  pre- 
sented principally  in  the  form  of  periodograms , precluding  direct  comparison 
with  the  time  series  generated  by  the  model.  Variance  computations  for  these 
results,  however,  suggest  a sigma  of  about  six  dB  which  is  in  reasonable  agree- 
ment with  model  results.  In  addition,  it  is  qualitatively  obvious  that  the 
periodograms  for  the  intermediate  range  hydrophone  fall  off  more  rapidly  with 
frequency  than  does  the  long  range  periodogram,  supporting  the  prior  suggestion 
that  propagation  fluctuations  osci  1 late  more  rapidly  with  increasing  range. 

Quantitatively,  these  periodograms  show  a strong  semi-diurnal  (12. 48 
hr)  spectral  component  (hypothesized  to  arise  from  tidal  effects)  which  is,  of 
course,  not  present  in  the  model.  At  the  other  extreme,  they  cut  off  at  10 
cycles  per  hour,  a limit  imposed  by  the  rate  (once  per  three  minutes)  at  which 
the  original  data  was  sampled.  Most  of  the  power  in  the  spectrum  appears  to 
lie  well  below  the  10  cph  cut-off.  This  result  appears  to  contradict  model 
predictions;  a one  minute  relaxation  time  corresponds  to  a frequency  bandwidth 
of  about  60  cph  which  Is  much  larger  than  the  less  than  0 -ph  bandwidth  suggested 
by  the  data.  However,  firm  comparison  of  these  time  scales  requires  subjecting 
the  model  results  to  the  same  data  reduction  algorithms  (filtering,  sampling) 
as  the  data,  and  this  has  not  been  done  here. 

This  point  may  be  emphasized  by  turning  to  the  Eleuthera  experimental 
results  where  a somewhat  more  precise  comparison  with  model  results  is  possible. 
Here  we  have  available  a reasonably  good  sample  time  series,  which  is  reproduced 
in  Figure  4.43.  To  compare  this  time  series  with  the  model  predictions,  level 
crossing  rates  (number  of  crossings  per  hour)  were  estimated  for  levels  in 
5 dB  steps,  from  5 dB  down  from  the  time  series  maximum,  to  30  dB  down.  The 
results  are  shown  In  Figure  4.44. 

Such  a plot  compares  both  the  rate  of  fluctuations  observed  (on  the 
ordinate)  and  their  magnitude  (on  the  abscissa).  We  observe  that  the  model 
generated  level  crossing  rates  are  within  30  percent  of  those  obtained  experi- 
mentally (approximately  19  crossings  per  hour  versus  approximately  14  crossings 
per  hour)  while  the  widths  of  the  two  curves  are  substantially  the  same.  We 
conclude  that,  In  this  case,  the  correspondence  between  model  predictions  and 
experimental  results  is  rather  good. 


1 

I 
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FIGURE  4.43  ACOUSTIC  FLUCTUATIONS-ELEUTHERA  EXPERIMENT 


4.4  SIGNIFICANCE  AND  CONCLUSION 

Aside  from  specific  details  of  the  simulation  results  presented  above, 
several  general  features  are  notable.  We  may  summarize  these  as  follows: 

§ The  rate  at  which  propagation  fluctuations  occur,  for 
a fixed  source  and  receiver,  increases  (the  relaxation 
time  decreases)  as  the  acoustic  frequency  increases 

• There  is  some  evidence  that  fluctuation  rates  increase 
with  increasing  source- receiver  separation 

• Both  predicted  and  observed  fluctuations  tend  to  a 
maximum  value  of  about  5 to  6 dB,  a value  which  is 
characteristic  of  phase  random  multi -path  propagation. 

One  set  of  runs  shows  evidence  that  this  value  decreases 
with  range  at  short  (less  than  40  nmi)  separations  but 
the  current  runs  are  insufficient  to  determine  whether 
this  feature  changes  with  source-receiver  depth,  frequency, 
or  environment 

• Computed  relaxation  times  range  from  one  to  two  minutes 
up  to  about  30  minutes.  Values  at  220  Hz,  at  ranges 
of  100  to  200  nmi,  are  about  five  minutes  in  both  the 
Bermuda  and  Atlantic  environments,  and  on  the  order  of 
one  to  two  minutes  In  the  two  400  Hz  cases  examined. 

All  of  these  points  are  of  potential  significance  to  the  performance 
of  passive  narrowband  ASW  sonars.  The  magnitude  and  rate  of  received  signal 
and  noise  are  critical  factors  in  determining  the  performance  of  such  systems. 
While  propagation  effects  are  certainly  not  the  only  contribution  to  these 
parameters,  these  effects  can  be  important,  under  proper  conditions,  to  the 
point  of  dominating  system  performance. 

Conventionally,  sonar  fluctuation  parameters  are  taken  to  be  indepen- 
dent of  range  and  essentially  independent  of  frequency  over  the  operating 
range  of  any  particular  sonar.  Empirical  fits  to  low  frequency. long  range, 
sonar  performance  suggest  relaxation  times,  for  the  detection  process,  of 
hours  rather  than  minutes.  The  possibility,  suggested  by  the  results  above, 
that  fluctuation  parameters  may  vary  with  both  frequency  and  range  could  be 
of  considerable  importance  in  optimizing  sonar  design  and  use.  Crudely 
speaking,  the  cited  propagation  fluctuation  characteristics  imply  relatively 
better  performance  with  Increasing  range  and  frequency,  than  is  predicted 
by  conventional  calculations.  For  example,  to  the  extent  that  such  effects 
dominate,  sonars  designed  to  reach  a given  range  could  have  lower  gains,  and 
operate  at  higher  frequencies,  then  would  otherwise  be  necessary. 

Of  course,  many  phenomena  besides  those  treated  here  can  produce  or 
modify  observed  fluctuations  on  any  particular  sonar.  There  are  for  example 
time  fluctuations  Introduced  when  source  or  receiver  move  through  spatial 
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propagation  structure.*  There  are  fluctuations  in  ambient  noise  which  arise 
from  an  ensemble  of  sources  (hence  from  a superposition  of  many  spatial  ana, 
perhaps,  frequency  regimes)  whose  properties  may  differ  significantly  from 
the  fluctuation  characteristics  of  a single  source.  Additional  fluctuations 
arise  from  variations  in  the  properties  of  the  sources  themselves.  Other 
mechanisms,  not  treated  here,  may  also  perturb  acoustic  propagation. 

Finally,  the  sonar  system  Itself  modifies  the  properties  of  an 
Incoming  signal  so  that,  ultimately,  the  observed  fluctuations  depend  on 
the  properties  of  the  system.  Long  integration  times  will  remove,  or  smooth, 
relatively  rapid  fluctuations.  High  resolution  systems  can  introduce 
fluctuations  by  responding  to  changes  In  parameters  which  would  not  be 
resolvable  on  less  powerful  systems.  This  is  not  to  suggest  that  the  trends 
observed  in  this  study  are  not  of  interest,  but  only  that  realistic  assess- 
ment of  sonar  performance  must  integrate  such  results  Into  a more  comprehensive 
sonar  model.  As  such  models  are  developed,  estimates  of  transmission  statistics 
will  play  a central  role. 


* When  the  spatial  structure  arises  from  multipath  Interference,  such 
fluctuations  are  closely  related  to  those  discussed  above.  That  Is,  they 
also  arise  from  varying  phase  relations  among  multipaths.  However,  in  this 
case  the  time  scale  Is  set  by  platform  speed  rather  than  rate  of  change 
of  medium  properties. 
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V.  SUMMARY  AND  CONCLUSIONS 


Understanding  and  predicting  the  performance  of  low  frequency  passive 
sonars  requires  being  able  to  describe  not  only  mean  values  of  parameters  that 
determine  system  performance,  but  statistical  variations  as  well.  Of  these 
parameters,  propagation  loss  has  a central  role,  since  it  helps  determine 
the  properties  of  both  the  observed  signal  field,  and  the  observed  noise 
field  at  a receiver.  Of  a number  of  phenomena  which  can  generate  propagation 
fluctuations,  this  report  has  concentrated,  in  particular,  on  how  internal 
wave  phenomena  affect  acoustic  propagation  between  a stationary  source  and 
receiver. 


A model  has  been  developed  which  simulates  multipath  propagation 
through  an  internal  wave  field.  The  model  depends  on  building  up  sample 
functions  of  an  internal  wave  field  from  a specified  spectral  description, 
and  propagating  acoustic  energy  through  this  field  using  ray  trace  techniques. 
The  internal  wave  calculation  uses  a power  spectrum  proposed  by  Garrett  and 
Munk,  and  assumes  an  isotropic  distribution  of  internal  wave  propagation 
directions.  While  the  Internal  wave  spectrum  used  is  a single  universal 
function,  the  acoustic  propagation  calculation  is  based  on  propagation  para- 
meters which  are  generic  to  particular  environments  of  interest;  hence  the 
sensitivity  of  propagation  loss  under  different  environmental  conditions  to 
internal  wave  effects  may  be  tested. 

Qualitative  features  of  model  results  match  typical  experimental 
observations  of  CW  transmission.  Of  these  the  most  striking  are  frequent 
deep  fades  In  the  propagation  amplitude,  and  the  character  of  related  phase 
fluctuations  which  take  place  relatively  slowly,  but  with  large  (many  cycles) 
amplitude.  In  two  detailed  comparisions  of  model  results  with  data, one  shows 
very  good  agreement  In  both  the  magnitude  and  rate  of  fluctuations.  An 
apparant  anomoly  in  fluctuation  rates  observed  in  the  other  comparision  may 
result  from  differences  in  processing  and  display  between  simulated  results 
and  experimental  data;  this  case  requires  further  examination. 


' 
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Model  results  suggest  a number  of  trends  which  may  be  significant 
from  the  standpoint  of  assesing  sonar  performance.  These  include: 

• An  apparant  increase  in  fluctuation  rates  with 
increasing  acoustic  frequency 

• A possible  increase  in  fluctuation  rates  with 
increasing  source-receiver  seperation 

• Standard  deviations,  for  propagation  fluctuations, 
of  5 to  6 dB  under  all  conditions  except  short 
(approximately  40  miles)  source-receiver  seperations 

• No  obvious  strong  dependencies  on  environment,  or 
on  source-receiver  depth,  in  the  cases  tested. 


All  these 
sonar  performance, 
the  characteristics 
tuations  In  critica 
predictions  in  the 
model  uses  simulati 
additional  simulati 
from  model  results 


trends  are  potentially  important  in  the  assessment  of 
but  their  significance  is  hard  to  assess,  not  least  because 
of  any  sonar  system  may  themselves  affect  signal  fluc- 
1 ways.  Thus  an  important  next  step  is  to  evaluate  model 
light  of  specific  sonar  processing  concepts.  Because  the 
on  techniques  this  can  be  done  in  two  ways:  directly  through 
ons,  or  indirectly  by  using  statistical  parameters  derived 
in  sonar  system  analyses. 


A possible  application  of  considerable  interest  would  be  to  evaluate 
limits  imposed,  by  internal  wave  activity,  on  Inter  Array  Processing  (IAP). 
This  technique  depends  on  having  coherent  signal  arrivals  at  two  widely  sepa- 
rated sensors.  This  coherence  can  be  degraded  or  destroyed  by  perturbations 
Introduced  by  Internal  waves.  Adaptations  of  the  model  developed  under  this 
study  could  readily  simulate  this  effect.  Other  possible  applications  include 
driving  particular  signal  processing  algorithms  with  signal  time  dependencies 
derived  from  the  model,  and  assessing  limitations  on  sonar  array  beamforming 
capabilities  caused  by  coherence  degredation  across  a wavefront.  In  each 
case  the  model,  as  developed  In  this  study,  offers  substantial  potential  for 
investigating  limits  of  sonar  performance. 
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APPENDIX  A 

MATHEMATICAL  FORMALISM 


BASIC  PROPERTIES 

In  this  section  we  outline  some  of  the  mathematical  formalism  which 
underlies  the  description  of  internal  wave  activity  given  in  Chapter  IV.  The 
treatment  Is  based  largely  on  the  1974  paper  of  Garrett  and  Munk.  5 In  par- 
ticular we  use  the  Garrett-Munk  estimate  of  the  power  spectrum  of  internal 
waves.  We  also  adapt  the  Garret-Munk  solutions  for  the  normal  modes  of 
internal  waves.  These  are  based.pn.the  assumption  that  the  internal  wave 
bouyancy  frequency  (the  "Brunt- Vai sal  a"  frequency)  varies  exponentially  with 
depth.*  Estimates  of  the  effect  of  Internal  wave  activity  on  acoustic  propa- 
gation are  based  on  ray  trace  estimates  of  unperturbed  propagation  parameters, 
and  an  extension  of  an  approximation  introduced  by  Porter  Spindel  and  Jaffee.6 

We  treat  internal  waves  in  the  ocean  as  a superposition  of  traveling 
waves  of  the  form: 


6(a,  w) 


iLjaifl  art.  »)  e1'01 


(1) 


This  is  equivalent  to  assuming  (see  Appendix  B,  Equation  12)  that  the  depth 
derivative  of  the  logarlthlm  of  sound  speed  varies  exponentially  with  depth. 
This  characterization  is  meant  only  to  capture  general  trends,  and  real 
sound  speed  profiles  frequently  depart  from  it  radically.  While  we  use 
the  Garrett-Munk  internal  wave  model  functions  derived  from  this  approx- 
imation, acoustic  ray  trace  and  transmission  parameters  were  based  on 
estimates  of  actual  sound  speed  profiles,  at  sites  of  interest. 
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where: 


<5  Is  vertical  displacement  due  to  Internal  waves 
a Is  a two  dimensional  (horizontal)  wave  vector 
to  Is  Internal  wave  frequency 

Z Is  the  vertical  wave  function  associated  with  mode  (cT,to) 

A Is  the  complex  amplitude  (a  random  variable)  of  mode  (a,<o) 

nQ  Is  the  peak  Internal  wave  bouyancy  frequency. 

The  total  displacement  at  a given  location  r,  at  time  t Is,  formally, 
just  a sum  over  the  constituent  modes.  That  Is: 


6(r,  2,  t) 
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_ The  mean  square  of  A Is  the  mean  energy  carried  by  mode  (a,<o).  If 
E (a,<o)  is  defined  as  the  power  spectrum  (mean  energy  density)  associated 
with  Internal  waves  then: 


{|A|2}  * - - E (a,w)  6to 

2ir|a| 

where  62  a 6co  Is  the  size  of  a cell  containing  one  mode.  In  wave  number  fre- 
quency space.  Garrett  and  Munk  propose,  for  E (ot.to),  the  form 
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where  v * 20  rr  /w2  _ ^ 

E = 2it  x 10'* 

b * mean  bouyancy  frequency  depth 
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nQ  3 peak  bouyancy  frequency 

u>j  3 Inertial  frequency  3 (2  c^jes)  sin  ( 1 at ) 3 .000145  sin  (lat) 

a 3 v/  ( b n ) 
max  ' o' 

(In  this  formula  all  frequencies  are  circular  (e.g.,  radians/sec)  and  the 
scaling  factors  used  by  Garrett  and  Munk  are  displayed  explicitly.) 

Garrett  and  Munk  found  that  the  discrete  modal  structure  of  Equation 
(2)  predicts  spatial  coherencies  for  Internal  wave  fields  which  are  larger 
than  those  observed  experimentally.  They  suggest  that  this  reflects  a more 
complex  modal  distribution  and  propose  approximating  this  by  replacing  the 
discrete  structure  In  (2)  with  a continum  In  (a,cj)  space.  In  this  case  the 
summation  In  (2)  is  replaced  by  an  integral  over  this  continuous  distribution. 

We  can  retain  the  discrete  structure  of  (2)  by  partitioning  the  Integral 
Into  bins  within  which  a and  u>  are  nearly  constant.  In  each  bin  the  phase 
factor  Is  nearly  constant,  so  that  the  Integral  Is  essentially  over  the 
terms  a A.  This  Integral  Is  again  a Gaussian  random  variable  with  mean 
square  equal  to  the  sum  of  the  mean  square  values  of  Z and  A separately. 

This  leads  to 

6(?,z,t)  3 2 £££  A'(a,w)  e7’^  ' 7 * wt)  (4) 

o,w  no 

where  {Z2}*5  Is  the  root  mean  square  value  of  Z over  an  (cT.oj)  group,  and  the 
sum  now  Is  over  averaged  groups  rather  than  single  modes.  A’  is  again  a 
Gaussian  Random  amplitude,  with  mean  square  energy 

{A'2}  3 E(oT,co)  <S2a'  5to 1 (5) 

where  62a'  6w'  Is  the  size  of  the  cell  over  which  averaging  takes  place. 

The  solutions  for  single  modes, Z(a,u>)  as  developed  by  Garrett  and 
Munk,  exhibit  a strongly  oscillatory  depth  dependence  which  varies  rapidly 
from  mode  to  mode,  and  the  solutions  are  sensitive  to  detailed  features  of 
the  depth  dependence  of  oceanographic  variables.  Garrett  and  Munk  point 
out  that  the  mean  square  average  of  Z Is  a much  smoother  and  less  sensitive 
function  of  the  form 


uj.|  < a)  < n(z)  (6) 
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where 
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To  connect  the  vertical  excursions,  6,  to  acoustic  propagation,  we 
observe  that  these  oscillations  are  mixing  up  water  of  different  temperatures. 
This  In  turn  affects  the  speed  of  sound.  In  Appendix  B it  is  shown  that  the 
resultant  change  In  sound  speed  at  depth  z may  be  written  as 

yc 

4C(«)  • - (-g£)  n2(«)  e (7) 

where  g Is  the  acceleration  of  gravity,  u is  a proportionality  constant,  and 
c01s  the  unperturbed  local  sound  speed.  The  change  in  sound  speed  is  equivi- 
lent  to  a change  in  the  local  Index  of  refraction  given  by 

N,(r,.,t)  -4§-  - - (-£)  n*(i)  « (8) 

a c0  9 


Inertial  frequency  = 2 sin  (latitude) 

Brunt  - VSisSlS  frequency  = {-  Ij}15 

acceleration  of  gravity 
density  of  seawater 
peak  value  of  n( 2) 


While  perturbations  In  the  refractive  index  can  occur  throughout  the 
water  column,  the  Garret-Munk  formulation  suggests  that  the  major  contributions 
tend  to  He  In  a well  localized  region  at  the  base  of  the  thermocline.  Porter, 
et.  alstake  advantage  of  this  by  Introducing  an  approximation  which  treats  the 
sum  of  all  the  variations  in  Index  throughout  the  water  column  as  taking  place 
In  a thin  layer  at  the  base  of  the  thermocline.  That  is,  they  assume  that  at 
the  base  of  the  thermocline  there  Is  a thin  slab  whose  total  acoustic  thickness 
differs  from  that  of  the  unperturbed  water  column  by 

/Bottom 

Na(r,z.)  dz  (9) 
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In  this  context  It  Is  straightforward  to  treat  the  effect  of  internal 
waves  on  a ray  trace  model  of  acoustic  propagation.  Rays  which  do  not  cross 
this  thin  Internal  wave  layer  are  unaffected.  For  those  rays  that  do  cross, 
reflections  from  the  layer  are  Ignored  (because  the  change  In  index  of  re- 
fraction Is  small),  any  horizontal  displacement  of  the  ray  Is  Ignored  (because 
the  slab  Is  thin)  and  any  change  In  direction  Is  Ignored  (because  the  boundaries 
of  the  slab  are  parallel).  The  remaining  effect  Is  a perturbation  in  phase 
as  the  ray  crosses  the  perturbed  region,  and  this  is  given  by 
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(10) 


a*  * Hr  Na(rj’  V/,cos  ejl  (10) 

where  X Is  the  acoustic  wavelength  of  the  ray,  ej i s the  angle  (from  the 

vertical)  at  which  the  ray  crosses  the  slab,  and  r . and  t,  are  the  location 
and  time  of  the  crossing.  J J 

In  an  actual  propagation  situation  there  can  be  several  rays  con- 
necting source  and  receiver,  and  each  ray  can  cross  the  Internal  wave  layer 
more  than  once.  Each  time  a ray  crosses  It  receives  an  additional  phase 
perturbation.  Let  the  positions  where  ray  k crosses  the  internal  wave  layer 
depth  be  rj,  and  let  the  corresponding  crossing  times  be  t - t*  (t  the  arrival 

time  at  the  receiver).  Then  the  total  phase  perturbation  for  ray  k is: 


¥ S 

A J 
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To  connect  this  to  acoustic  propagation,  note  that  the  total  acoustic 
transmission  from  source  to  receiver  Is  given  by 


I (If.  t) 


£ Ek 
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where  and  are  the  underturbed  amplitude  and  phase,  respectively,  of  ray 

k,  and  <Sdk  Is  given  by  (11).  Note  that  the  perturbing  phase,  as  given  by  (11) 
Is  explicitly  a function  of  time,  and  implicitly  a function  of  source  range 
(since  this  determines  how  many  times  each  ray  crosses  the  layer). 


EXPLICIT  REPRESENTATIONS 

Equation  12  Is  the  formal  result  required,  but  we  are  not  yet  in 
posltloji  to  use  It.  In  order  to  evaluate  (ll)we  need  an  explicit  expression 
for  N»(r,t).  This  requires  evaluating  the  Integral  In  (9)  (including  the  2 
dependence  In  6),  which.  In  turn,  requires  supplying  an  explicit  representation 
forn(z).  The  first  step  simply  involves  combining  (8),  (9)  and  (4),  giving: 
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where 


Bottom 


f dz  n2(z)  (Z2(z)  J*5 


As  noted  above  the  Garrett-Munk  formalism  assumes  that  n(z)  is 
exponential  In  form,  that  Is: 

n(z)  * n0e"Z/b 

Adapting  this  assumption,  and  using  the  expression  for  {Z2(z)}  given  by 
(6 ) we  have 


(cj2-  u^.2)1* 
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where  the  upper  limit  of  the  Integral  in  (15)  is  set  by  the  range  over  which 
(6)  is  non-zero. 

Combining  (18)  with  (14),  and  using  this  result  in  (11),  yields  the 
following  expression  for  the  phase  change  along  path  k: 
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Part  of  th£  summations  In  (19)  may  be  evaluated  in  closed  form.  The 
propagation  vector  a Is  a two-dimensional  vector  in  the  horizontal  plane.  If 
we  take  the  X-axis^  of_our  co-ordinate  frame  to  be  the  propagation  direction 
of  Interest  then  a • r = axx  and  (19)  becomes 
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where  A"(ax,  u>)  » £ A‘(a,  “>)  (21) 

ay 

A"  Is  again  a complex  mean-zero  Gaussian  Random  variable.  In  Appendix 
C we  show  that  Its  variance  Is  given  by: 
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where 
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The  crux  of  using  (12)  to  evaluate  the  acoustic  Intensity  (I)  Is  to 
find  suitable  ways  of  evaluating  (20).  Quantities  such  as  Ek>  «Sk,  r!j,  t!j,  0j, 

can  be  obtained  readily  from  a suitable  ray  trace  pro£ram.  The  crucial  problems 
are  to  select  a suitable  mesh  size  for  the  sum  over  (a.u),  and  to  handle  the 
random  variables  A"(a,w). 


To  take  the  second  problem  first,  our  technique  Is  to  Implement  a 
simulation  approach  In  which  values  for  the  real  and  Imaginary  part  of  A" 
at  each  (a,w)  mesh  point,  are  selected  randomly  from  a mean  zero  Gaussian 
distribution  with  variance  given  by  (22).  For  each  set  of  random  numbers 
so  selected,  (20)  and  (12)  can  be  evaluated  at  all  values  of  IT  and  t which 
are  of  interest.  The  result  is  one  realization  of  the  random  function  I (r,t). 
Iterating  this  procedure  can  build  up  an  ensemble  of  such  functions  from  which 
relevant  statistics  (mean,  variance,  autocorrelation,  distribution  function, 
etc.)  as  functions  of  R and  t,  can  be  generated. 


The  question  of  selecting  a mesh  size  for  ax  and  w Is  somewhat  more 

subtle.  We  approach  the  subject  as  follows:  If  a uniform  mesh  Is  chosen  (e.g., 

6ax  and  5co  constant)  then  for  any  given  value  of  j,  equation  (20)  Is  a double 

Fourier  series.  That  Is,  ignoring  the  sum  over  j,  we  have  a function  which 

Is  periodic  In  range  with  period  2ir/|6aJ,  and  periodic  in  time  with  period 

2ir/6w.  This  periodicity  Is  quite  clearly  an  artifact  generated  by  the  way 

the  terms  In  equation  (2)  were  grouped  to  form  equation  (4),  rather  than  a 

real  effect.  To  keep  It  from  affecting  the  results  It  is  necessary  that  the 

period  be  longer  than  all  ranges  and  times  of  Interest.  That  Is,  If  we 

select  some  maximum  range,  R_.v,  and  some  maximum  time  Tm  . of  Interest,  we 
. max  ma  x 

require 
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The  prescription  for  carrying  out  the  calculation  Is  then  as  follows 

a)  Pick  a mesh  according  to  (3) 

b)  Randomly  choose  a value  for  the  real  and  Imaginary  parts 
at  A"  at  each  mesh  point 

c)  Carry  out  the  sum  In  equation  (20),  and  evaluate  equation 
(1?.)  for  each  IT  and  t of  Interest 

d)  If  desired  repeat  b)  and  c)  to  build  up  statistics. 
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APPENDIX  B 

SOUND  SPEED  - BOUYANCY  FREQUENCY  RELATIONSHIPS 


in 

BASIC  RELATIONSHIPS 

We  wish  to  ascertain  how  vertical  mixing  of  the  water  column  affects 
the  local  speed  of  sound.  The  following  treatment  Is  based  on  a discussion 
given  by  Munk  and  Zacharlaseni 3 We  observe  that  sound  speed  (c)  depends  on 
temperature  (T),  salinity  (S),  and  pressure  (P)  so  that 

I 3c  „ 3T  + 8 li  3P 

C 32  32  6 32  Y 32 

(1) 

where  a,  B,  and  3 are  the  fractional  partial  derivatives 
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The  temperature  gradient  may  be  divided  Into  two  parts;  an 
representing  the  thermal  gradient  that  would  be  present  In 
and  the  remaining  ("potential")  gradlant.  That  Is 

adiabatic  componet  .. 
an  adiabatic  ocean. 

3T  m 3Ta  + 3TP 
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so  we  can  write 
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where 
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Is  the  fractional  sound  gradient  expected  In  an  adiabatic  isohallne  ocean; 
It  Is  approximately  the  gradient  observed  In  the  deep  ocean. 


When  an  element  of  the  water  column  Is  displaced  from  its  equilibrium 
depth  by  mixing,  it  carries  with  It  Its  original  value  of  salinity  and  Tp,  but 
takes  on  the  local  values  for  P and  Ta.  Thus  an  element  which  has  been  dis- 
placed a distance  6 has  associated  with  It  a sound  speed  which  differs  from 
the  unperturbed  speed  at  that  depth  by: 


5c(2)  „ _1 
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C 32 
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The  gradient  of  C_  depends  in  an  approximately  linear  fashion  on  the 
gradients  of  Tp  and  salinity.  Similarly  the  Internal  wave  bouyancy  frequency 
depends  on  the  same  partial  derivatives  through  a relationship  of  the  form 
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where  Ty,  the  "Turner  Number"  Is  defined  as 
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In  terms  of  Tu  there  Is  a relationship  between  ^ and  n*(2) 


which  may  be  expressed  as 


B-2 


(9) 


1 

c 


. a 
g 


nJ(a) 


where 

a 1 + (a8/ab)  T 

u ■ - — y~t t E 

v 


(10) 


Thus,  combining  (6)  and  (9)  we  have 

Y’  * _ g 6 (ID 


which  expresses  the  fractional  change  In  sound  velocity  In  terms  of  the 
bouyancy  frequency  n(a). 

Finally,  we  note  that,  from  (4)  and  (9),  we  have 
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EVALUATION  OF  CONSTANTS 

The  proportionality  factor  p,  which  appears  above.  Is  expected  to 
vary  with  both  depth  and  location.  Typical  values  quoted  In  the  literature 
vary  from  11.6  to  24.5;  for  the  calculations  reported  here,  a value  of  18 
was  used. 


In  equation  A16  the  approximation 
n(a)  * nQ  e~z/b 

Is  Introduced.  From  equation  (2)  we  have,  also,  that 

"(*)  * i H'f  <13> 

This  relation  may  be  used  to  fit  for  nQ  and  b,  using  the  local  velocity 
profile  at  any  location  of  Interest  as  an  Input.  In  the  work  reported  here, 
this  was  done  by  fitting  to  the  sound  speed  slope  at  the  base  of  the  thermo- 
cllne  / 3c  | _ 3c  i \ and  the  depth  (h)  of  the  sound  channel  axis 
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Resulting  values  of  b are  on  the  order  of  1 to  2 Km;  of  n , about  .005  Rad/Sec. 
Specific  values  for  the  cases  treated  In  Chapter  IV  are  given  In  Table  4.1. 


APPENDIX  C 

SUMMATION  OVER  TRANSVERSE  WAVE  NUMBER 


Equation  (21)  reads 
A"  (a,w)  * ^ A'  (a.u) 


where  oi  Is  a two-dimensional  wave  number  vector  with  cartesian  components 
(ax»av),  A'  (a,w)  Is  a complex  mean  zero  Gaussian  random  variate  with  variance 
gl9enyby  (Equation  A5). 

{A12}  * -L—  E (a.w)  62  o'  6u'  (2) 

2ir|a| 

and  the  sum  Is  over  a set  of  closely  spaced  grid  points  In  a,u  space.  The 
values  of  A'  at  the  various  grid  points  are  uncorrelated  so  that 


r (ax,oa)2}  - J){A'  (a,c 


T!  1 E(a,oj)  6ax'  6ay'  6u> 


ay  2tt  ( ot  ( 


From  Equation  (A3),  E(a,co)  has  the  form 


E(a,to) 


E(u)  e\,1n(u)  1 |a|  < » M 


otherwise 


Thus,  E(a,oj)  does  not  depend  explicitly  on  a and  we  may  write 


E(<o)  6a'  6co 
{A"  (a  ,<o)2}  xr* 
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^ pi 


where  the  limits  on  the  sum  in  (5)  are  determined  by  the  range  of  |a|  as 
given  In  (4).  Using 

H a ’'a  2 + a 2 we  have 

A J 

min  . i | . max  u 

ay  i I ®y I 1 ay  « 


where  a. 
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The  range  of  ay  covers  both  positive  and  negative  values.  However 
the  sum  In  (5)  Is  symmetric  so  that  we  may  write 


max 

E(w)  6a  1 5aj  ay 
(A”  <v  „)«)  • 2^ £ 


a ' ■ a™1" 
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Finally  we  take  advantage  of  the  smallness  of  the  mesh  spacing 
( 6oty ' ) to  replace  the  sum  In  (6)  with  the  equlvllent  Integral,  so  that 

amax 

{A"  (a.w)2}  = ^ 6tJ  ■■  /am1n  day  H ay*  + <*x2 
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